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, , Abstract 

The interpretation of some recent measurements of spin-dependent voltage for which the electric 

conduction does not play a role rises some new fundamental questions about the effects of spin- 
CD 

dependent heat currents. A two spin-channel model is proposed in order to describe the effect of 
out-of-equilibrium spin-dependent heat carriers in electric conductors and insulators. It is shown 
that thermal spin-accumulation can be generated by the heat currents only over an arbitrarily long 
distance for both electric conductors or electric insulators. The diffusion equations for thermal spin- 

o 

O accumulation are derived in both cases, and the principle of its detection based on Spin-Nernst 

^ effect is described. 
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The role played by thermoelectric effects in the context of giant magnetoresistance (GMR) 
has long been recognized [IrE] and exploited for possible applications in spintronics. Indeed, 
spin-dependent Seebeck and spin-dependent Peltier measurements are direct generalizations 
of giant magnetoresistance experiments performed in the presence of temperature gradi- 
ent, that can be described in the framework of the two spin-channel model for conduction 
electrons [TUTU] . Starting from the coupled transport equations of thermoelectricity, the in- 
troduction of the spin-dependent conduction coefficients a± and the spin-dependent Seebeck 
cross-coefficients S$ suffices to describe most of the spin-dependent thermoelectric effects in 
electric conductors [TTHTS] . However, recent experiments pioneered by the group of Uchida 
et al. [T9H21] and still not well understood [TlSTfTl] opened the way to a new class of spin- 
tronics effects in which no electric current is flowing and only spin-dependent heat currents 
are present in the sample. Such experiments are performed in electric conductors in the 
so-called non-local geometry (Fig. 1) at a long distance from the current injection, or in 
electrical insulators (Fig. 2). 

The goal of this work is to propose a phenomenological description of spin-accumulation 
produced by spin-dependent heat currents. The model first describes the spin-thermocouple 
effect that allows spin-accumulation to be generated by a temperature gradient in conduc- 
tors over arbitrarily large distances. This spin-dependent version of the thermocouple effect 
(that allows a temperature difference to be measured with a voltmeter in an open circuit), is 
corollary of the Seebeck effect (that allows an electric current to be generated with a temper- 
ature difference in a closed circuit). Both effects are deduced from the same thermoelectric 
transport equations with different boundary conditions and different power dissipation. The 
model is then applied to the case of electric insulators within the two spin-channel model 
of heat conduction. Finally, the transverse spin-dependent thermocouple effect occurring 
in the spin- Hall or spin-Nernst detector is also described. In all cases, the corresponding 
diffusion equations for the spin-accumulation are derived. 

For a non-isothermal ferromagnetic material - that can be an electric conductor or an 
electric insulator - the system under interest is composed of a statistical ensemble of out- 
of-equilibrium heat carriers of two different species: the heat carriers that carry a up spin 
■f and the heat carriers that carry a down spin J,. In the case of electric conductors, the 
spin-dependent heat carriers are also electric carriers, while for electric insulators, the spin- 
dependent heat carriers are spin-dependent excitation modes (e.g. magnons). In both cases 
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(II) J e = 

FIG. 1: Schematic view of the spintronics non-local device with temperature gradient. Branch (II) 
under interest constitutes an open circuit from the point of view of the electric charges ( J e = 0) 
but not from the point of view of the heat currents. The voltage is measured between points A 
and B in branch (II), which contains ferromagnetic/non- ferromagnetic interfaces, or between the 
points C and D of the lateral electrode. 
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FIG. 2: Schematic view of non-local device with ferromagnetic insulator. Heat current is injected 
through the device and spin-accumulation can be produced in the lateral electrode. The potential 
due to spin-dependent Nernst effect is measured transversally along the y direction. 

however, the two sub-systems of t and j. heat carriers are described by the corresponding 
local chemical potential fi^ and fii, assuming the hypothesis of local equilibrium [25j. In 
line with the description of the giant magnetoresistance [7H9], we define the thermal spin- 
accumulation by the difference of the two chemical potentials Afi = ^ — fi± of the heat 
carriers. In the framework of the non-equilibrium thermodynamic theory, the bivaluated spin 
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variable ^ defines an internal degree of freedom JTUJ |25j |26j. The spin-relaxation occurring 
during the transport process is then treated as a chemical reaction that transforms heat 
carriers from the up channel to the down channel or inversely. The rate ip th of this reaction 
is a flux of spin variable in the spin space. Since this reaction changes the density of the heat 
carriers of the spin channels, the corresponding chemical affinity is A/i. The Onsager relation 
that links the flux ip th to the force A/i takes the same form as for electric spin-accumulation 

[na mill]: 

4 th = L th Afi, (1) 

where we have introduced the phenomenological Onsager coefficient L t h that describes the 
spin-dependent relaxation. At stationary regime (i.e. for the state of minimum entropy 
production of the system), the internal energy densities and u_i obey the following con- 
servation equations: 

^ = -V.4-^ + ^/2 = 

^ = -V.Jj 7 + ip th + 11/2 = 0, 
where 1Z describes the dissipation out of the system (e.g. by radiation) and ip th describes 
the spin relaxation from one heat conduction channel to the other. 



In the case of electric conductors and in a first approximation, the heat carriers and the 
electric carriers are the same, so that the chemical potentials of the heat conduction channels 
is also that of the electric conduction channels. The well-known thermoelectric transport 
equations relate the electric currents J| and the heat currents to the corresponding forces 
V/^ and VT [161 HU [271 : 

je = _ a l^ h _ a Ay Ti (3) 

4 = X t VT - 5* (4) 

where is the spin-dependent Fourier coefficient and 11^ = cr^S^T is the spin- dependent 
Peltier coefficient. The electrochemical potential reads jl^ = [i c h,x + e ^ where e is the 
charge of the electric carriers and V the electric voltage. We assume that the temperature 
dependence of the transport coefficients is negligible. 

The condition for the partial equilibrium Jj = is then given by 

= -S X VT (5) 
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Eqs. ([5]) can be reformulated with introducing the total Seebeck coefficient So = S[ + S±: 

V (/2t + = -So VT (6) 

and the Seebeck asymmetry coefficient AS = S^- — 64.: 

V(A^) = -A5VT. (7) 

Eq. ([6]) is nothing but the equation of the thermocouple. Equation ^ is, in contrast, a direct 
consequence of the presence of spin-dependent Seebeck coefficients, and describes the spin- 
dependent thermocouple effect. Accordingly, it is expected that a gradient of temperature 
produces a spin- accumulation gradient without electric current and without non-local spin 
injection. The usual thermocouple effect allows the spin-accumulation generated by the 
temperature gradient to be measured: inserting Eq. §f§ in to Eq. (|6| we obtain: 

V(/i t + ^) = ^V(A yU ). (8) 

Eq. (|8p is able to account for the observed Spin- Seebeck effects measured without electric 
current (Jf = Jf — 0) at a distance arbitrarily large in the configuration of Fig. 1. The 
integration between the points A and B leads to the spin-dependent potential measured: 

Vs-V^r^^- (») 



l A AS dx 

On the other hand, Eqs. (J3]) to Eq. (J8j) can be generalized with the introduction of ten- 
sorial transport coefficients {cr±}ij, {$±}ij (hj = {x, y, z) ) where the non-diagonal elements 
i 7^ j account respectively for the Hall effect and the corresponding Nernst effect in ho- 
mogeneous isotropic materials [25J. The detection of the transverse potential difference Vdc 
performed on the right electrode in Fig. 1 and Fig. 2 is described by the transverse spin- 
dependent thermocouple effect, or spin-dependent Nernst effect. If AS xy 7^ 0, the projection 
of Eq. ([7]) over the Oy axis gives, after integration: 

v °- v ° = L Da w dv= £ As ^ dv - (io) 

Eq. Q and Eq. §LQ§ are analogous to that of the usual spin-accumulation calculation 
for the giant magnetoresistance effect, with the difference that the spin- accumulation is pro- 
duced by the heat current instead of the electric current, and that the spin-accumulation 
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A/i(x) and A/j,(y) are no longer governed by the usual diffusion equation of the spin- 
accumulation generated by electric current injection [7HIU]. 

Defining the total heat flux along the x axis (see Fig. 1 and Fig. 2) Jq = Jl + J^ and he 
spin- dependent heat flux A J q = J* — JJ, we have from Eq. (El): 



(11) 



U u oi 2e \ ox ' ox / 



where /i = /Sf+A|,j A = A t +A 4 , AA = A t -A 4 , n = n t +Il 4 , and (3 = (Il t - ILJ / (Il t + . 

Note that for Jq = 0, the flow A J 9 is a "pure spin- current' [T2l [13] that cannot be 
distinguished from its electric homologue AJ e = Jf — Jf with Jq = Jf + J| = (and more 
generally, a "pure spin-current" is independent of the nature of the spin carriers involved). 

Inserting Eq. (11) and Eq. Q into Eq. we deduce the diffusion equation for A/i: 



^"If -^fJ?^ 2 ' (12) 
where we have introduced the thermal spin-diffusion length for conducting materials 



and the constants d = 2e (AA - A /3) / (n (l - /3 2 )), C 2 = 2f31le/ (n (l - /3 2 )). 

Assuming a uniform temperature gradient and negligible external dissipation 1Z ~ 0, 



the solutions of Eq. (12) have the form A/i(x) = A/x(0) (ae x ^ lth + be~ x ^ lth ^j where a and b 
are constants imposed by the boundary conditions. It is the same solutions as that of the 
usual spin-accumulation responsible for GMR effects, with the difference that the diffusion 
length is defined with the help of the Peltier asymmetry coefficient (5 = AIl/Ilo instead 
of the asymmetry of the electric conductivity, and L t h plays the role of the usual spin-flip 
Onsager coefficient [101 [11] . 

In the case of spintronic devices with electric insulators (Fig. 2), there are no electric 
carriers and the chemical potentials /i^ are that of the corresponding heat carriers only 
(magnon, spin-dependent phonon, etc). The heat current is composed of a drift part Jj ri f t $ = 
A^VT and a diffusion part JdiffX = — At^/z^ 1 , where is the corresponding Onsager 
coefficient related to the diffusion coefficient D$ = kTC^/n th , where n th is the density of the 
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heat carriers (for the sake of simplicity, we assume isotropic diffusion). In the configuration 
of Fig. 2, the transport equations reads: 



ji — \ QT _ r a ^\ h 
J t ~~ A t dx M dx 

j q = \ ®L — C d ^ 

X v dx 4- dx 



(14) 



The diffusion equation for the chemical potential difference (i.e. spin accumulation) is de- 
duced by inserting Eq. (14) and Eq. Q into Eqs. (j2]): 



d 2 Afx th Afi th ~ d 2 T ~ 



dx 2 l 2 n dx 2 
where we have introduced the thermal diffusion length for electric insulators 



(15) 



hi 



\ 



Co 1 - &< 



4L 



tii 



(16) 
2(AA - 



where C = + C± and (3 = — C±) / Cq. On the other hand, C\ 
XoP)/ (C (l - (3 2 )) and C 2 = 2^/((l - (3 2 )C ). 

This spin-accumulation in the electrical insulator is transferred to a conducting electrodes 
through the condition of continuity of the heat currents at the interface. In that case also, 
the final profile of the spin-accumulation throughout the insulator /conductor interface is 
similar to that of the usual GMR spin-accumulation under electric current. 



In conclusion, the two spin-channel model has been applied to spin-dependent heat trans- 
port, in a circuit that is closed from the point of view of heat currents, and open for the point 
of view of electric currents. The result of the model shows that the heat current is able to 
play the same role as the electric current in a usual "local" spin-valve system, whatever the 
spin-dependent heat carriers considered. In particular, spin-injection and spin-accumulation 
can be performed with spin-dependent heat currents only at an interface. The model ex- 
plains the recent observations of spin-dependent voltage measured in electric conductors of 
electric insulators, in the configurations depicted in Fig. 1 and Fig. 2. 
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